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Abstract 

Numerical simulations are performed to study the finite temperature phase tran- 
sition in the SU(2) Higgs model on the lattice. The strength of the first order phase 
transition is investigated by determining the latent heat and the interface tension on 
L t = 2 lattices. The values of the Higgs boson mass presently chosen are below 50 
GeV. Our results are in qualitative agreement with two-loop resummed perturbation 
theory. 



1 Introduction 

The standard calculational method for the study of the symmetry restoring electroweak phase 
transition [|TJ is resummed perturbation theory [Q, |3|, Q . The highest known order in the gauge- 
and scalar quartic coupling is g A , A 2 , combined with a high temperature expansion ||. In the 
Higgs phase with broken symmetry perturbation theory is expected to work well for not very 
high Higgs boson masses, since the couplings are small. The perturbative calculations can be 
extended to the vicinity of the symmetric phase for A <C g 2 . In the high temperature symmetric 

*On leave from Institute for Theoretical Physics, Eotvos University, Budapest, Hungary. 



1 



phase, however, the situation is similar to high temperature QCD. Therefore, irreparable in- 
frared singularities occur which prevent a quantitative control of graph resummation |J . Since 
the calculation of physical quantities characterizing the phase transition requires the knowledge 
of both phases, there is a priory no reason why perturbation theory could provide a quantitative 
treatment of the electroweak phase transition for physically relevant couplings. 

An obvious non-perturbative method for the study of the symmetry restoring phase tran- 
sition is numerical simulation on the lattice. Omitting, for simplicity, fermions and the U(l) 
gauge field one is left with the SU(2) Higgs model describing the interaction of a four-component 
Higgs scalar field with the SU(2) gauge field. After pioneering works ]7|, |8[], recent numerical 
simulations contributed a lot to the understanding of the finite temperature behaviour of the 
SU(2) Higgs model M. Another line of research has also been initiated recently [ID, II]. In this 



approach one studies the three-dimensional effective theory obtained by dimensional reduction 
in the high-temperature limit. 

An important issue is the order and strength of the phase transition. For instance, the pos- 
sibility of electroweak baryogenesis [12j requires a strong enough first order electroweak phase 
transition [|13] (for a recent review see [14|]). Up to now the four dimensional numerical simula- 
tions H have been concentrated on relatively heavy Higgs boson masses, where the properties 
of the phase transition are technically difficult to disentangle. In the numerical approach based 
on dimensional reduction indications have been found that the phase transition is stronger 
than given by perturbation theory |I3| . In the present letter we investigate this question in 
the four dimensional SU(2) Higgs model for lighter Higgs masses, below 50 GeV. Since this 
region of parameters of the minimal standard model is already excluded by experiments, our 
present scope is merely theoretical because we would like to check the validity of some other 
theoretical approximation schemes, e.g. resummed perturbation theory. We plan to extend this 
investigation to heavier Higgs boson masses in future papers. 

The present letter is our first short communication. Technical details will not be included 



here, but will be postponed to a subsequent publication with more numerical data |16 



2 Physical parameters 



The lattice action of the SU(2) Higgs model is conventionally written as 

S[UM=^{l-\TrU^ 



P i 



+ E ^Tr fa+p.) + A 



(1) 



Here U x ^ denotes the SU(2) gauge link variable and tp x is a complex 2 <g> 2 matrix in isospin 
space describing the Higgs scalar field and satisfying 



r 2 ^r 2 



(2) 



The bare parameters in the action are (3 = A/g 2 for the gauge coupling, A for the scalar 
quartic coupling and k for the scalar hopping parameter related to the bare mass square /Zq by 
H% = (1-2A)/k-8. 

In order to fix the physical parameters in a numerical simulation one has to define and 
compute some suitable renormalized quantities at zero temperature. The renormalized gauge 
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coupling can be determined from the static potential of an external SU(2) charge pair, measured 
by Wilson loops. The physical Higgs mass Mh can be extracted from correlation functions of 
quantities as the site variable 

R x = -Tr [<p+(p x ) = p 2 x , (3) 
or, using (p x = p x a x , the link variables 



L v = \ Tr (ft+pUx^x) , L a ee ^Tr (a+ +p U xll a x ) . (4) 



The W-boson mass Mw can be obtained similarly from the composite link fields (r, k = 1,2, 3) 

W rk ee ^Tr (r r a+ + -U xk a x ) . (5) 

Since we are interested in the study of the symmetry restoring phase transition on lattices 
with temporal extension L t = 2, we have to determine the renormalized parameters at critical 
points for the L t = 2 lattices. In order to have a renormalized gauge coupling near the physical 
value g\ ~ 0.5 we choose (3 = 8 frFfl . As stated before, we would like to have lighter Higgs boson 



masses than studied in ||. Therefore we have chosen the two values A = 0.0001 (referred to as 
low) and A = 0.0005 (referred to as high). In these cases the critical hopping parameters turned 
out to be near k 1 ° w = 0.1283 and K^ l9h = 0.1289, respectively. For the numerical simulation of 
the zero temperature (T = 0) system we took 12 3 • 24 lattices, which turned out to be large 
enough for the present L t = 2 simulations. (Runs on 16 3 ■ 32 are planned in near future fT6| .) 



The numerical Monte Carlo simulations have been performed on the Quadrics Q16 ma- 
chines at DESY. In order to decrease autocorrelations we applied a mixture of heatbath and 
overrelaxation algorithms with an optimized ratio, which depends on lattice sizes and bare 
parameters. An important improvement could be achieved by an overrelaxation algorithm in 



the Higgs field length (p) variable which has been developed and tested by two of us |18| . The 
heatbath algorithm for the gauge field can be taken over from pure gauge theory ||19|| . For the 
Higgs field one can start a heatbath algorithm at our small A values by generating an optimized 
Gaussian part of the single spin distribution and then correct by an accept-reject step for the 
remaining quartic piece | 20fl . Away from the critical k such combination of updating sweeps, 



which we call "complete sweep", almost entirely removes the autocorrelation. On our 12 3 • 24 
lattices the integrated autocorrelations for different characteristic quantities are between 2 and 
20 complete sweeps. The smallest integrated autocorrelations are shown by Wilson loops, the 
largest ones by variables proportional to the length of the Higgs field. In the metastability 
region with tunneling between the two phases we still have the "supercritical slowing down", 
however, in the simulations presented here this does not play a role. In the two-coupling simu- 
lations of sections 3 and 4 the observed integrated autocorrelations were maximally of the order 
of 100 to 300 complete sweeps. 

Performing 160000 complete sweeps on 12 3 • 24 lattices and using an exponential + constant 
fit for the timeslice correlations of the quantities in (0-^) one obtains the following masses in 
lattice units: 

aM l ° w = 0.236(7) , aM^ i9h = 0.262(9) , 

aM$» = 1.059(24) , aM^ gh = 0.427(8) . (6) 

Here, as everywhere in this paper, statistical errors in last digits are given in parentheses. 
The estimates of the statistical errors always come from dividing our simulation data into 
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independent subsets and performing the fits in these subsets in order to obtain the variance of 
the fit parameters. This procedure is particularly simple to implement on a Quadrics parallel 
computer, where in most cases several independent lattices are simulated. In our case we 
started Liischer's random number generator |2T| differently on each node, leading to statistically 



independent samples on the independent lattices. 

In general we do not yet attempt to estimate our systematical errors, because they are 
presumably dominated by the present limitation to L t = 2 lattices for the high temperature 
simulations. In the continuum limit one has to consider L t ^> 1, which is possible, but can be 
demanding [ |i~6f . 

For the mass ratio Rhw = Mh/Mw = V8Xr/9r (H) gives 

R l ™y = 0.222(12) , Rx§ = 0.614(32) . (7) 

With My/ ~ 80 GeV this means that the low point is at Mu ~ 18 GeV and the high point at 
Mh — 49 GeV. On the L t = 2 lattices we have aT c = 0.5, therefore (|6]) also implies 

T l ° w = 0.472(11) • M w , Tf 9h = 1.171(22) • M w . (8) 

The static potential has been extracted from r ® t Wilson loops with 1 < r < 6 and 
1 < t < 12 on 12 3 • 24. We determined every such Wilson loop after transforming the gauge 
configuration to temporal gauge. The t-dependence has been fitted by three exponentials in 
order to determine the large-t asymptotics. The potential can be well fitted by a Yukawa- 
potential with lattice corrections, as discussed in jr?j. Taking a form 



V(r) = -je~ Mr + C + D- G(M, r) , (9) 

with G(M, r) being the difference of lattice Yukawa-potential and continuum Yukawa potential, 
we obtained with \ 2 of order one: 

A low = 0.0336(3) , M low = 1.112(15) , C low = 0.066453(6) , D low = 0.0343(23) , 

A Mgh = o.03434(7) , M hi9h = 0.4272(21) , C high = 0.090768(18) , D high = 0.0352(8) . 

' (10) 

The excellent agreement of the masses M in ( ]T0| ) with Mw in (0) and the nearly equality 
of the fit parameters A and D show that the static potential is well understood, and can be 
used for the definition of the renormalized gauge coupling. One could take for a definition 
g R = 16nA/3, but we prefer here not to use a global fit. Instead, we adopt for the definition of 
the renormalized gauge coupling a procedure similar to the one recently proposed in pure gauge 



theory P^| . Namely, starting from eq. (35) in ref. [IT] we define an r-dependent renormalized 
gauge coupling, which we then interpolate to a distance r = M w x . The result is, this time 
including also an estimate of the systematic errors as the last entry in parentheses: 

g^ low = 0.5476(5 + 75 + 10) , g 2 R Mgh = 0.5781(4 + 88 + 3) . (11) 

The first two entries are the statistical errors from the Wilson loops and M\y, respectively. As 
one can see, the errors are dominated by the statistical errors of Mw used as input in this 
analysis. This can be easily improved by better statistics, if later on required. Otherwise the 
errors are quite small, and the renormalization compared to the bare value g 2 = 0.5 is moderate. 
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3 Latent heat 



The pressure (P) is continuous, therefore the latent heat (i. e. the discontinuity of the energy 
density Ae) can be calculated from the discontinuity of the quantity 5 = e/3 — P. Since 5 is 
given by [|3| 



where r = \og(aMw) , and besides L v in (|J) we used 




g. = (p£-l) a , P^l-^TrfV (13) 



The vacuum contribution is not subtracted in Q12D , but in the latent heat it cancels and we 
obtain 

The partial derivatives are taken here along the lines of constant physics (LCP's), which 
can be defined by keeping constant the value of the mass ratio Rrw an d the renormalized 
gauge coupling g R . For weak bare couplings one can estimate the change of g 2 = 4//3 and 
Ao = A/(4k 2 ) by integrating the one- loop renormalization group equations 

dg\r) 1 



dr 16tt 2 



~g A + O(Xl,Xlg 2 ,X g\g e 



d\ (r) 1 



96A ° + h 9 * ~ 9X092 + ° (A °' A °^' A °^ 4 ' ^ 



(15) 



dr 167T 2 

Starting from our low or high critical point of the L t = 2 lattice one can obtain the corresponding 
L t = 3,4, 5, . . . points by integrating flTBD from r = to r = log(3/2), log(4/2), log(5/2), . . .. 
For instance, the L t = 3 critical point corresponding to our low point (J3 = 8, A = 0.0001, k = 
0.12830) comes out to be at {(3 = 8.147, A = 0.000111) and the one corresponding to the high 
point (p = 8, A = 0.0005, k = 0.12887) at (j3 = 8.147, A = 0.000507). 

Of course, we still need another equation for the change of the critical k = k c . The critical 
hopping parameter can be determined by different methods. We found it convenient to apply 
a two-coupling method: on a long lattice in one spatial direction (say, z-direction) the lattice 
is subdivided into two equal halves with two different hopping parameters in such a way that 
the first half with K\ is in the symmetric phase, the second half with K2 in the Higgs phase. 
This means K\ < k c < k^. To be sure that this situation is stable we require that both halves 
stay in their phases at least for 20 autocorrelation times. Going to larger lattices k 12 can be 
chosen closer to k c at the order of inverse lattice volume. From the 2 • 16 2 • 128 lattice at the 
low point and 2 ■ 32 2 ■ 256 lattice at the high point the best estimates are, respectively, 

<°Z=2 = 0-12830(5) , «££=2 = 0.12887(1) . (16) 

These are rather well reproduced, together with the volume dependence on smaller lattices, 
by the critical k obtained from the one loop gauge invariant effective potential of the length 
square of the scalar field V e //(p 2 ) 0, ^5[ |TB|. Therefore one can use V e //(p 2 ) for an estimate 



of K c ,L t =3,4,5,- -- By a quadratic interpolation of the L t = 2,3,4 critical points one obtains the 
following estimates for the derivatives of k along the LCP's: 

8k 1oW Q^-high 

-—- = -0.00064(6) , — — = -0.00111(13) . (17) 
or or 



5 



The discontinuities in ([14]) can be determined at the low point on large enough lattices 
directly at k — k c , because the metastability is very strong and the configurations stay for a 
very long time in the phase they started from. At the high point a linear extrapolation can be 
used from the intervals k < K\ and k > k 2 , where no tunnelings to the other phase occur. On 
our 2 • 16 2 • 128 and 2 • 32 2 • 256 lattices at the low and high points, respectively, the finite volume 
corrections of order (L x L y L z )~ 1 are already smaller than the indicated statistical errors: 



AL io W 



-21.84(4) , AQ l < 



low 



high 



AL hi 9 h = _ L 026(12) , AQ 
Combining this with (|T^)-(|l5D we arrive at 




7^7 = 1-68(17) 



-643(2) 



A plow 



0.010557(15) 



-8.02(13) , APu i9h = 0.000629(8) . (18) 



> A \high 

jrA =0.125(19) 



(19) 



These results have been obtained on L t = 2 lattices alone. In the continuum limit L t — > oo, 
therefore later on one has to improve them by simulations on finer lattices with L t > 2 [|T6 
For a first qualitative comparison with two-loop resummed perturbation theory H see fig. 1 



4 Interface tension 

Besides the latent heat another important physical quantity characterizing the strength of the 
phase transition is the interface tension a between the two phases, which can be determined on 
the lattice in several different ways. In the present letter we concentrate on the two-coupling 



method proposed by Potvin and Rebbi | 26fl , where an interface pair is enforced between the 
two halves of the periodic lattice. Since we have three bare parameters, in principle one can 
choose any of them (or some combination of them) to be different in the two halves. As in 
the latent heat the contribution of the </)-link L v dominates, we decided to split the value of 
k multiplying L v in the action. On an L t ■ L x ■ L y ■ L z lattice with L z ^> L x ^ yt , and in our 
case L t = 2, L x = L y , the two halves in the ^-direction have hopping parameters K\p which 
are slightly below and above the critical one: Ki < k c < k 2 . The interface tension between the 
states with k x and k 2 is given in lattice units by 

3 <7(Ki, K 2 ) = L z \ r dKL$>(K, K 2 ) - f 2 dKL®(K U k)\ , (20) 

where L^' 2 \k, k') denote the expectation values of the yj-link in the two halves if the hopping 
parameters are k and k', respectively. 

In order to obtain the physically interesting interface tension one has, of course, to perform 
the limits L XiV ^ z — > oo and An = — > 0. For a given lattice extension, however, 
An cannot be arbitrarily small, because if the difference in free energy density becomes small 
tunneling into the other phase can occur and the interfaces disappear. The presence of the 
interfaces can, however, be monitored to ensure the applicability of (^0|). For small An the 
integrals in ( |20"D can be well approximated by the average of the values of the integrand at the 
two end points. Besides, for equal arguments we obviously have L^(k,k) = L@\k,k). This 
gives 

a 3 a ~ L z Ak (L$\ Kl , k 3 ) - I$\k u «i) + L^(k 2 , k 2 ) - L<?\ Kl , « 2 )) . (21) 
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In the limit An —>■ the contributions with equal hopping parameters in the two halves do not 
contribute. The non-zero contributions can be determined by the expansion 

l£\ki, k 2 ) = — h h + ai(«i - k c ) + 0(ki - k c ) 2 , 

« 2 ) = ~° 2 + b 2 + a 2 (K 2 - K C ) + 0(k 2 - n c f . (22) 

This formula replaces (|20|) in the region where ft 2 ) depends linearly on «i i2 . For Ak — > 
a finite volume estimator of the interface tension is 

a 3 a = L z ( Cl + c 2 ) . (23) 

The advantage compared to fl2"0|) is that no numerical evaluation of integrals is necessary. 

In our numerical simulations several lattice sizes up to 2 • 16 2 • 128 at the low point and up to 
2 • 32 2 ■ 256 at the high point were exploited. In fact, the very long extensions in the z-direction 
have been chosen just for the purpose of also determining the interface tension. The measured 
L^' 2 ^ values were fitted in a carefully chosen /t-interval around the critical points in ( |T6"D by the 
three-parameter forms in ( p2]) . At the low point we used the intervals 0.1279 < K\ < 0.1282 and 
0.1284 < k 2 < 0.1287 with seven points each. At the high point we used 0.12880 < m < 0.12886 
and 0.12888 < k 2 < 0.12894 with seven points each. The number of complete sweeps for 
measurements with an optimized ratio of heatbath to overrelaxation was between 20000 and 
40000 per point. The obtained result from the fits in different sub intervals is: 

(» \ low / " \ high 

=0.84(16), i^j = 0.008(2). (24) 

Here aT c = | has been used in order to convert from lattice to physical units. About half of the 
indicated errors come from statistics. The other half is a crude estimate of systematic errors, 
which is due to the choice of subintervals for k. In general, choosing a subinterval closer to k c 
gives smaller statistical errors and also somewhat smaller ci )2 . The quoted results are obtained 
from the more distant subintervals. 

This dependence on the chosen fit interval is a kind of rounding-off effect which is partly 
due to the interaction of the two interfaces. A similar dependence of ct(ki, k 2 ) on k 12 appears 
in (|20| ) if «i )2 are too close to k c . This ^-dependence is expected to be substantially reduced for 
L z —>■ oo when the two interfaces do not interact. Of course, the transversal extensions L XjV also 
have to be taken to infinity. In our case comparisons to smaller lattices show in both points 
that within the given statistical errors L x = L y are large enough. In principle, if L xy z are 
larger one can consider smaller An values, but this is beyond the scope of the present analysis. 

We also tried several other methods for the determination of a. The results are in general 
qualitatively similar. A more detailed evaluation of the different methods will be included in 
a subsequent publication [jT6fl . A comparison of (|24D with two-loop resummed perturbation 
theory || is shown in fig. 2. 
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5 Discussion 



The outcome of our numerical studies on L t = 2 lattices is a qualitative agreement between two- 
loop resummed perturbation theory and four dimensional lattice simulations in the investigated 
range of Higgs boson masses, below 50 GeV. Note that at Mh — 18 GeV the agreement between 
our data and perturbation theory becomes even better if, instead of T c , M H is used to set the 
scale. Both latent heat and interface tension show that the first order symmetry restoring 
phase transition becomes substantially weaker for increasing Higgs boson mass (see eqs. ( |T9| ) 
and (|24"D ). As it has been emphasized throughout this paper, although the L t = 2 results may 
very well be qualitatively right, for a quantitative determination of the properties of the phase 
transition further numerical simulations on L t > 2 lattices are necessary. 
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Figure captions 

Fig. 1. Comparison of the results of Monte Carlo simulations with two-loop resummed 
perturbation theory ||. The latent heat Ae is shown as a function of the Higgs boson mass. 
The two curves correspond to the two slightly different gauge couplings in ([11]): the full line to 
the low value, the dashed line to the high one. 

Fig. 2. The same as fig. 1 for the interface tension a . 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9405021v2 
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